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Abstract 

The method of A-operators developed by S. Derkachov, G. Korchemsky, A. Ma- 
nashov is apphed to a derivation of eigenfunctions for the open Toda chain. The 
Sklyanin measure is reproduced using diagram technique developed for these A- 
operators. The properties of the A-operators are studied. This approach to the 
open Toda chain eigenfunctions reproduces Gauss-Givental representation for these 
eigenfunctions. 

1 Introduction 

This work was inspired by the article [DKM] devoted to the Separation of Variables (SoV) 
method for XXX-mode\. The main idea of this method is to find an integral transfor- 
mation such that eigenfunctions of quantum integrals of motion in new variables becomes 
the product of functions of one variable [Skj . If everyone of these functions satisfies the 
Baxter equation, then the initial multivariable function becomes an eigenfunction. The 
kernel of this transform is called a transition function and can be constructed as consec- 
utive appUcation of operators Ak{u): Aj^{ji)An_i{'J2) ■ ■ ■Ai(jn). Every operator A^^u) 
is an integral transformation, which maps a function of A; — 1 variables onto function of 
k variables. The properties of the transition function can be translated to the properties 
of these operators (see section [5l). 

The transition function for the A^-particle periodic Toda chain was obtained in the 
works [Gutzj , [Sk] , |PGj [KLlj . The transition function in this case is proportional to the 
eigenfunction of the (A^ — l)-particle open Toda chain, with a factor depending on the 
A^-th particle coordinate. We apply methods of the paper [DKM] to obtain these eigen- 
functions as a product of A-operators reproducing the Gauss-Givental formula |GKLOj . 
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This form of eigenf unctions of the open Toda chain leads to an integral represen- 
tation that appeared first in [Givj employing a different approach. Recently it was 
interpreted from a group-theoretical point of view using the Gauss decomposition of 
GL{N,'R) [GKLUj . Therefore, this integral representation of the eigenfunctions for the 
open Toda chain is called a Gauss-Givental representation. 

The method of a triangulation of the Lax matrix described in [PG] was used in |DKM] . 
We also use a triangulation, which is implemented by a gauge transformation parame- 
terized by variables yo,...,yj\f. In the periodic case one has to impose the condition 
Vo = Un and the method described in |PGj produces Baxter's Q-operators |Baxj for the 
periodic Toda chain model. Following [DKMJ we impose a different boundary condition: 
yo — oo, i/AT +00 to construct A-operator. Thus A-operator and Baxter's Q-operator 
for the periodic Toda chain correspond to the different choice of the boundary conditions 
in the method of triangulation of the Lax matrices. 

To describe the construction of eigenfunctions for open Toda chain we develop a kind 
of the Feynman diagram technique similar to one exploited in |DKM] . It allows to reduce 
calculations with kernels of A-operators to simple manipulations with diagrams. 

The article is organized as follows. In section [5] we recall a definition of the open Toda 
chain model following |KL2t [KL3] . Section [H] is devoted to a description of eigenfunctions 
in terms of the product of A-operators and formulation of a diagram technique developed 
m |DKMj . In section H we use this technique in order to prove that eigenfunctions satisfy 
an orthogonality condition. As a by-product of this calculation we obtain a Sklyanin 
measure, which is necessary to prove a completeness condition. Section is devoted to 
algebraic properties of A-operators. 

2 Open Toda chain model 

The quantum A^-particle open Toda chain is a one- dimensional model with the exponential 
interaction between the nearest particles. The hamiltonian of the system is equal to 

^=2E^^' + E^""^ (2.1) 

n=l n=l 

d 

where pn = —ih— — is an operator of momentum for the n-th particle. Due to the 

OXn 

translational invariance, the total momentum 

N 

P = J2pn (2.2) 

n=l 

commutes with the hamiltonian, i.e. it is also an integral of motion. There are func- 
tionally independent integrals of motion for this system. It is relevant to use the i?-matrix 
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formalism to find them. First of all, introduce the Lax operator of the Toda chain 

Lniu)=(''-J: n = l,...,iV, (2.3) 



^ — e^'" 

and monodromy matrix for the A^-particle Toda chain 

where « is a spectral parameter. 

The following recurrent relations, which are direct consequence of this definition, will 
be useful below: 

An{u) = {u~ pn)An^i{u) + e-"^C^_i(ii), (2.5) 

CNiu) = -e"^A,v-i(«), (2.6) 

Ar,{u) = iu~pN)AN^,iu) - e^--i-^-^^_2(n). (2.7) 

These relations show that Ajv(m) and Cn{u) are polynomials in u of degree and — 1 
respectively. Analogously, Bj^{u) and Dn{u) have degree A^ — 1 and N — 2. 

The monodromy matrix satisfies to the quantum i?TT-relation 

R{u - v) {Tn{u) ® /)(/ ® Tn{v)) = (J ® Tn{v)){Tn{u) ® /) R{u - v) (2.8) 
with the rational i?-matrix 

R{u) = I^I+-V, (2.9) 

u 

where P is a permutation matrix: Vij^ki = ^u^jk- 

Rewriting fl2.8p by entries one obtains, in particular, the relation 

An{u)An{v) = An{v)A^{u). (2.10) 

This means that A]\f{u) is a generation function of the integrals of motion of the 
integrable system with A^ degree of freedom. Explicit calculations show that these are 
integrals for open Toda chain model: 

N 

An{u) = ^(-l)^^-^i/,, (2.11) 

Ho = l, Hi = P, H2 = ^P^-H, (2.12) 

[Hk,H,] = 0. (2.13) 

By virtue of (12.131) there exist common eigenfunctions of the integrals correspond- 
ing to the eigenvalues Ek- They are defining by the equation 

AN{u)ijE{x) = ajv(n; E)tPEix), (2.14) 
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where 

N 
k=0 

Eq = 1, E = {El, . . . , En), X = {xi, . . . , xn)- Representing eigenvalues Ek as symmetric 
combinations 

Ek= Yl ^n---l^. (2-15) 

j\<...<ik 

of real variables 7 = (71, ... , 7Ar), one can rewrite equation (12.141) as follows 

N 

An{u)iP^{x) = Y[{u - -fj)4j-y{x). (2.16) 

We also require the solutions of the equation (I2.16P to be rapidly decaying in region 
Xi ^ Xj+i. This asymptotic properties of solutions correspond to the condition that the 
i-th particle is situated on the left of the {i + l)-th particle. 

3 Eigenfunctions of the open Toda chain 

In this section we shall find eigenfunctions of the open Toda chain defined in the previous 
section by the equation (I2.16p . This equation is equivalent to the system of equations 

A^(7i)V'7(^) = 0, j = l,...,N. (3.1) 

The eigenvalues (12.151) are invariant under the permutations of 71, . . . ,7Ar. Therefore, it 
is reasonably to require the invariance of eigenfunction under these permutations, which 
we shall call the Weyl invariance: 

^o-(7)(a^) = i^-yix), for all a G Sn, (3.2) 

where Sn is a permutation group and (7(7) = (7cr(i), • • • ,7o-(Af))- 

It is sufficiently to find a Weyl invariant solution of the unique equation 

AN{-fi)iJ^{x) = 0, (3.3) 

which will be a solution for the whole system (13.11) due to its Weyl invariance. 

To solve the last equation we shall consider a gauge transformation of the Lax operators 
K{u) = MnLn{u)M-\, n = l,...,N (3.4) 

by the matrices 

Mn=(^^y^ J), n = 0,...,N. (3.5) 
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The deformed A^-particle monodromy matrix is 

Tn{u) ^ ^^[^M = LM ■ ■■L,{u) = M^T^{u)M^\ (3.6) 

\Cn{u) Dn{u) j 

In particular, we have 

L„{u)2i = ie^'^iu -pn- «e^"-i"^" + ie^""^"), (3.7) 
Cn{u) = tey^'ANiu) + ey^+y°BN{u) + Cn{u) - le^^'DNiu). (3.8) 

Here Ln{u)2i is a lower off-diagonal entry of the matrix Ln{u). 
Let us consider the auxiliary equation 

Lniu)2lWniu) = 0, (3.9) 

which has the following solution 

Wn{u) = exp|^M(x„ - ?/„_i) - le^'^-i-^'' - -^e^""^"}- (3.10) 
It is clear that the function 

Wuix;y) = Ylwn{u) = exp ^|^m(x„ - y„_i) - -e^-i"^" - -e^""^"} (3.11) 

n=l n=l 

is a solution to the equation 

CNiu)Wuix;y)=0. (3.12) 
In the limit yo — ^ — oo, yjy +oo the formula (13.81) gives us the equality 



Aj^{u) = -i lim e-y^'CNiu). (3.13) 

Therefore, multiplying the equation fl3.12p by — ze~^^e^"''^°+^^\ taking the same limit as 
in (13.131) and setting m = 71 we arrive to the equation (13.31) with the solution ip-y{x) = 
= A^^(x; y), where 



Au{x;y)=Au{xi,...,XN;yi,...,yN-i)= lim e^"(^«+^^)l^„(x; y) 

N N-1 , Af-1 



exp<i -M(^a;„-^l/„) - -^^(e^" ^"+1 + e^" ^" 

n=l n=l n=l 



(3.14) 



Let An{u) be an operator with the kernel Au(xi, . . . , xn', yi, ■ ■ ■ , yN-i), i-e. 

{An{u) ■ f){x) = / (i?/A„(xi,...,a;7v;?/i,...,yAr_i)/(y). (3.15) 



Since |A„(x; y) \ decays at yi ±00 the integral in the right hand side of (I3.15P converges 
absolutely for the functions f{y) increasing sufficiently slowly in the infinity. Setting 
u = 7i in (I3.15P we obtain a solution to (13.30 . Since \Au{x;y)\ decays in the region 
Xi ^ Xj+i this solution has the requiring asymptotic properties. 
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Theorem 1. The following solution to the equation (\3.3h 

ip^^x) = (AAr(7i) ■ ■ ■ A2(7Af-i)Ai(7Af) ■ • • ■,xn), 



(3.16) 



where (Ai(7Ar) ■ = er^'^'^^^, is Weyl invariant, i.e. satisfies to the condition (13.21) . 

and, therefore, is a solution to the equation (12.161) rapidly decaying in the region Xi ^ Xj+i . 

Proof. It is sufficient to establish the invariance under the elementary permutations, i.e. 
we need to check the equality 



dy A^^_-^{Xi, . . . , XAr_„+2; Vl, ■ ■ ■ , yN-n+l)-A-f„{yi, ■ ■ ■ , UN-n+l] ^l, ■ ■ ■ , Zn-ti) 



-n+l 



dy A-y„iXi, . . . ,XN-n+2;yi, ■ ■ . ,2/7V-n+l)A^„_^(?/i, . . .^yN-n+l, Zi, 



ZN-n), 

(3.17) 



for n = 1, . . . , — 1. 





fig. Ic 



fig. Id 



Zk 



Xk+1 

fig. le 



To do it we shall use a diagram technique introduced in |DKMj . Let us denote 
the function I{xk,yk) = exp| — — e'^'=~^'= | by a line pictured in the fig. la, the function 

Ju{xk+i,yk) = exp|^u(a;fc+i - yk) - ^6^"-'"'+^^ by a fig. lb, the function Z^{xk) = 

= exp|— Mx^l by a fig. Ic, the function Z^^(xfc) = exp| — — -ux^j by a fig. Id and the 

i 

function YjyXk^i^ Zk) = ( 1 + e^^^^"^^ \ by a fig. le. 
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fig. 2 

The function fl3.16p can be represented in these graphical notations. For the case N = 5 
it is pictured in fig. 2, where bold bullets • signify that we integrate over corresponding 
variables. 

Lemma 1. The equalities represented in the figures 3a, 3b, 3c are valid. 




fig. 3a 





Proof. Integration over in the left hand side of fig. 3a yields 

o 

dykHxk, yk)J^„_,{xk+i,yk)I{yk, Zk).Lf„{yk, Zk^i) = et(Tn-i^fc+i-7n^.-i)x 



— oo 
+00 

—00 

'(7n-7n -l) 
p^fe _|_ p2fe-l \ 2h 



(3.18) 



X 



where K^(z) is a Macdonald function |BE2j . Interchanging jn-i and 7„ in 03.181] we 
obtain the expression for the integral in the right hand side of fig. 3a 



2gi(7na;fc+i-7„-i2fc_i) 



g — Xfc + l _j_ g— 2(; 



X 

(3.19) 



><^K7.-7n-i) 0V(e-'= + e-'=-i)(e-'=+i + e--o) • 

The ratio of f l3.18p and (13.191) is exactly equal to Y^^_^_^^{xk+i, Zk)Y~^^^_^^{xk, Zk-i). 

The equalities shown in the fig. 3b and 3c can be proved analogously. □ 

Let us continue the proof of the theorem [H It is shown in fig. 4. The left hand side 
of (13.171) after application of fig. 3c is refiected in this diagram. Then, using the fig. 3a, 
one can move the vertical line picturing the function Ky„_^_^„(xj, Zj_i) to the right, as 
shown in the figure. When this line has arrived to the right one can apply the fig. 3b. In 
each step the parameters 7^1-1 and 7„ are interchanged, and, eventually, one has the right 
hand side of flXTTj) . □ 
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UN-n+l 



fig. 4 



Substituting the expression f l3.15p for the operator An{'Jn) with the kernel (13.141) 
to (I3.16P one obtains the recurrent formula 



'07i,...,7]V (-'^l 5 ■ ■ ■ 1 -^n) 



dyi... dyN-i ^7i,...,7Ar^i (z/i, • • • , ?/7v~i) x 



X 



^ N N-l ^ N-1 

n=l n=l n=l 



(3.20) 



The consecutive applications of this formula allow to derive the following integral repre- 
sentation for the eigenfunctions of open Toda chain model 



^^(^ATl, . . .,Znn) 



i 2 



Af-1 n f . N N-1 n 

n n ] ^ [^^ 5Z + 5Z '^^^ ~ 5Z 

71=1 j=l ^ j=l n=l j=l 



_ ^^^g^Tij-2n-lj _|_ g^n-lj-^nj + l j I. (3-21) 

^ n=l j=l J 



This is a Gauss-Givental representation [GKLO] , |Giv] of the Toda chain transition func- 
tion. 



4 Integration measure 

In this section we shall check the orthogonality condition using the diagram technique. 
The normalization function which appears in this calculation coincides exactly with the 
Sklyanin integration measure using in the SoV method for the periodic Toda chain 
model [Sk], [KLl],[KL2], [KL3]. 



Theorem 2. The functions i^-yix) defined by the formula (I3.16p satisfy to the orthogonality 
condition 

j dxip^{x)ijY{x) = n'^{'y)SsYM{7,7')^ (4-1) 



where 



1 

m 



N 



is a symmetrized delta function and /i(7) is the Sklyanin measure 



-N 



7m - 7fc 

ih 



) 




n r 



\ ih / 



(4.2) 



k<m 



r/ie feot/i szc^es o/ equality (14.11) are understand as distributions with arguments 7^^, . . . , 7^, 
depending on the pairwise different parameters 71, . . . , 77V. 

Since the integrals of motion H^. are Hermitian operators, the set of the functions 
ip-yix) is complete [GVj . This means that any function f{x) belonging to the Hilbert 
space I/^(R^) can be represented as an integral 



for some summable function ^'(7). As a consequence we have the completeness condition 



Proof of the theorem [21 First of all, we need to obtain a diagram representation 
of i/j-yix) in order to calculate the integral in the left hand side of fl4.ip using the diagram 
technique. Let us to consider the diagram for iIj^{x), which are shown in fig. 2 for N = 5 
and in fig. 6a for = 3, and to implement the following steps. 

First step. The imaginary unit is contained only in the functions J and Z. To 
reduce the conjugation of whole function i/j^ix) to the conjugation of the functions Z 
we decompose Ju{xk+i,yk) into the product Z-^{yk)I{yk, Xk+i)Zu{xk+i) (fig. 5a). This 
corresponds to the transition from the fig. 6a to the fig. 6b. 

Second step. We replace all Zu{xk) by Z~^{xk) (fig. 5b) implementing the complex 
conjugation and arrive to the fig. 6c, in which the function i{j^{x) is pictured. 




(4.3) 




(4.4) 





fig. 5a 



fig. 5b 
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Third step. Since the functions Z~^{xk) are attached to only one point, namely to 
Xki one can turn these functions in the manner shown in fig. 5c and fig. 5d. It means that 
we can represent ■0-y(a;) by the fig. 6d. 

Fourth step. Now we replace the product Z-'^{xk)I{xk,yk)Zu{yk) by Ju{yk,Xk) 
(fig. 5e) to obtain fig. 6e. 




fig. 5c fig. 5d fig. 5e 



Fifth step. Finally, reflecting the fig. 6e with respect to a horizontal line we obtain 
the fig. 6f. 




fig. 6d fig. 6e fig. 6f 
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To obtain a graphical representation for the integral in (14. ip we should join the di- 
agrams shown in the figs. 6a for ip'^/{x) and 6f for i/j-yi^x) in the points Xi, . . . ,Xn and 
integrating over them. The diagram obtained like this in the case = 4 is pictured in 
fig. 8a. Further we shall simplify it using the following equalities. 

Lemma 2. 7/7' 7^ 7^ the equalities represented in the figures 7a and 7b are valid. 




fig. 7b 



Proof. The fig. 7a and 7b mean the following equalities, which are consequence of the 
integral representation of the Gamma-function 

+00 +00 

dxZ^^{x)I{x,y)J^'^{y\x) = j dx exp|^[(7fc - i^)x + 7^] - ^(e^"^ + e^^^') 

—00 



dxJ^^{x,y)I{y\x)Z^^{x)= / dx e^Y>[^[{lk - l'j)x - ikV] -^{e"" + e^' ^) 



e 



h 



ey + ef 



7fc - Iq 
ih 



h >ft r I — — ^ 

in 
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The poles of Gamma-function are non-positive integers. Since the variables 7^- and 7^ 
are not equal to each other and have only real values the right-hand sides of the equalities 
encoded in the figures 7a and 7b are defined correctly. □ 

The calculation of the integral in the left hand side of (14.11) can be reduced to a cal- 
culation of this integral in the domain 7^- 7^ 7^, j -|- A; 7^ + 1. In order to do it we need a 
smooth partition of unity which is symmetric over 7 with a component vanishing in some 
neihgbourhood of the plane 7^ = 7^, j + 7^ + 1. 

Since the variables 7^ are pairwise different one can choose a positive number e > 
such that 6e < miufc^; [7^ — 7;|. Let r]^{x) be a smooth function of one variable such that 
< ri^{x) < 1 for all x G M, r]^{x) = 1 for < e, //^(x) = for |x| > 3e. Let us set 

N 



7) = n {}~ ^'^^'j " ^''^ 



j,k=l 
j+k^N+1 

(^(7; 7) = - r, (4.5) 

where Sn is the permutation group and aj = (7cr(i), ■ ■ ■ ,1ct{n))- The formula (14.51) defines 
corectly a smooth function. To prove this it is sufficient to show that the denominator 
in the right hand side of this formula never vanishes. Since the functions 0(7'; (77) are 
non-negative we have to prove that for all 7' G there exists a permutation a G such 
that 0(7'; (J^) 7^ 0. Now fix 7' G and define this permutation a G Sn in the following 
way: let us consider a subset of subscript pairs ©(7') = {(j, k) \ — jkl < 3e}. It has 
the property following: if {j, k) G 6(7') and (j, /) G 6(7') then k = I. Indeed, if A; 7^ / we 
have a chain of inequalities 6e < |7a; — 7; | < |7a: — 7j- 1 + |7j- — 7; I < 3e + 3e = 6e, which leads 
to contradiction. This property means that there exist a permutation a G Sn such that 
i = a{k) for all (j, k) G ©(7')- Now let us show that (^(7'; a^) 7^ for a = a^^a, where 
a G Sn is the longest permutation: a{j) = N + 1 — j . Noting that k = a{N + 1 — j) for 
all {j, k) G ©(7') we see that [7^ — ^a{N+i-m) \ > 3e for all j 7^ m. Hence 

N N 

^(7'; ^7) = n n (1 - - 7.(^+i-m))) = 1 ^ 0. 

j=l m=l 

The function (14. 5 p has the following properties: 
1- Yli 0(7') ^1) = 1 (a partition of unity), 

2. if 17^ — 7fc| < e for some j, k such that j + A; 7^ + 1 then ^(7'; 7) = 0, 

3. if |7j- — 7Ar+i-j| < e for all j = 1, . . . , then 0(7'; 7) = 1. 

The property [T] is obvious. The property [2] follows from the fact that if [7^ — 7^! < e for 
some i, k such that j + A; 7^ A^ + 1 then (^(7'; 7) = 0. The property [2] follows from the same 
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fact. Indeed, fix 7' G and suppose that |7j- — "yN+i-j\ < e for all j = 1, . . . , A^. If a 7^ 1 
then there exist such j that a~'^{N + 1 — j) ^ N + 1 — j. Denoting k = a~'^{N + 1 — j) 
one derives {jj — 7cr(fe)| < and hence 0(7'; a^) = 0. In the other hand we have proved 
that there exists a E such that 0(7'; 0-7) 7^ 0, that is 0(7'; 7) 7^ 0. Hence 

0(Y; 7) 

Due to the property 1 and the symmetry of the sought integral the last one can be 
rewritten in the form 



where 



2^(7'; 7) = 0(7'; 7) / dxtp^{x)tpy{x) 



Therefore, since the Sklyanin measure is symmetric over 7 it is sufficient to derive 
the equality 

1 ^ 

ni';i) = f^-\i)j^U^h--iN^i-,). (4.6) 

The property [2] means that calculating X(7'; 7) one can use the lemma[2]for j + k 7^ + 1, 
because the poles of Gamma-functions will be nulled by function 0(7'; 7). The property [3] 
will be needed further. 

The expression X(7'; 7) is the integral represented graphically in the fig. 8a multiplied 
by 0(7'; 7). It can be calculated by induction. We integrate step by step over the boundary 
points connected with 71, j'-^. First, we integrate over the extremely left point and the 
extremely right point, i.e. over xi and xn, using the fig. 7a and 7b respectively, and we 
obtain the fig. 8b. Like in the fig. 4, the vertical line arising in the left side is moved 
to right, where it is annihilated by the line arising from the right side. This process 
exchanges 71 with 7^. After the integration the factor 

i[-i.\^n_i_i[) (4.7) 



ih J \ ih 

arises, and thus one obtains the fig. 8c, where 71, 7^ begin to be connected with another 
boundary points (boundary in sense of fig. 8a), Now we integrate over these points. In 
k-th integration (A; = 2, . . . , — 1) 71 exchanges with 7^ and j[ exchanges with 7^ and 
one has the factor 

rf^)rf^)rfl^)rf^). (4.8) 

\ in J \ in J \ in J \ in J 
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The A^-th integration leads to the factor 

(2^/1)5(7;- 7^) ■(2^^)5(7^v-7i) 



(4.9) 



After this the variables 71, 7J, 7Ar and 7^ disappear completely from diagram - they have 
gone away to the factors (14.71) . (14. 8 p and (14. 9p . The remaining diagram (the middle part 
of fig. 8e) is exactly the initial diagram, but for iV — 2 which is depend on 72, ... , 7Ar-i, 
72, ... , 7jv-i- Thereby we have obtained the recursion formula 



4>{i'^l) / dXi . . .dXN ^^^,...,-fj^{xi, . . . ,XN)i'^''y'{xi, . . . ,Xn) 



7i,-,7]vV 

*(y;7)r('^Vf^ix 



ih 



ih 



X 



N-l 
k=2 



7i -7fc\-pnfc -^AT^^fc -7l^^^l -7fc 



ih 



ih 



ih 



ih 



■(2^n)5(7i-7jv)x 



X {27Th)6{'y'^ - ji) / dx2 . . .dxN-i'ilJ^2,...,-yN-ii^2, ■ ■ ■ ,XN-i)'ipy2,...,'y'N_^{x2, ■ ■ ■ ,xn-i)- 



Continuing the calculation for the integral (14. ip in the same manner we obtain the fol- 
lowing result. 



0(7'; 7) / dXi . . .dXN i^'y^,...,-yN{Xl, . . . ,XN)ij^'^'{Xl, . . . ,Xn) 



Hi -,7) n 



k,j=i >- 

k+j<N 



ih 



ih 



N 



l[{27rh)5{i-^N+i-j). (4.10) 
-I j=i 





fig. 8b 



fig. 8c 



74V 
fig. 8e 



N 



The support of distribution Y[ ^{lj~lN+i-j) belongs to the domain I7' — 7Ar+i_j| < e. 



j = 1, . . . ,N, in which ^(7'; 7) = 1 due to the property El Substituting 7' = 7Ar+i_j in 
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the product of Gamma-functions in right hand side of (14.101) we derive (14.61) with the 
expression (14.21) for the Sklyanin measure. □ 



5 Properties of the A-operators 

Baxter's Q-operators Qn{u) described in |PGj for the periodic Toda chain satisfy the 
following properties. 

(a) These operators commute for different values of the spectral parameters: 
[Qn{u),Q^{v)] = Q. 

(b) They commute with the transfer matrix i^iu) = An(u) + Dm{u) of the periodic 
Toda chain model: [QN{u),tN{v)] = 0. 

(c) Q-operator satisfies the Baxter equation 

iN{u)QN{u) = i^QNiu + ih) + i'^Q^lu - ih). 

In this section the similar properties for the operators AAr(M) defined in section [3] will 
be established. 

Proposition 1. A-operator has the following properties: 
(l) ANiu)AN-iiv) = ANiv)AN-iiu), 
(ll) An{u)An{v) = {u- v)An{v)An_i{u), 

(hi) CAr(w)AAr(M) = i"^"^AAr(M - ik), 

(iv) Bn{u)An{u) = i^-'^AN{u + ih). 

Proof. The equality (l) was proved in the theorem [H Property (ll) have been 
proved in |GKLO] by the direct calculation. Below we give a simplified version of 
this property. Let us notice that (ll) is valid on arbitrary N — 1-particle eigenfunc- 
tion '07i,...,7jv-i(2/i; • • • ' Vn-i)- Indeed, in according with the equation (12.161) the left hand 
side of (ll) is 

An{u){An{v)^P^^^,„^^j^_^){xi, xn) = An{u)'iIj^^^,„^^^_^^^{xi, Xn) = 

N-l 

= {u-v)Yliu- -fk)^^u-,iN-iA^i^ ■■■,xn) (5.1) 

k=l 

while the right hand side is 

{u - v){AN{v)AN-i{u)'il)^^^,„^.f^_J{xi, . . . , Xn) = 

N-l 

= {u-v)Yliu- 7fc)(AAr(t;)'0^,,...,^^_J(a;i, ...,xn) = 
k=i 

N-l 

= {u-v)Yliu- 7fc)V'7i,-,7iV-i,^(^l; ■■■,Xn)- (5.2) 
k=l 
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Comparing (15. ip and fl5.2p one concludes that the equahty (ll) is vahd on the functions 
ip-yix). Due to the completeness of the system of these functions any function that belongs 
to the domain of definition of the operators An~i{u) and An{u) can be represented in the 
form (14. 3p . Therefore the property (ll) is valid on any function of — 1 variables where 
the action of the operators A]\i_i{u) and Ajv(f) is well defined. 

To prove (ill) we need the formula for an action of the operators Am{u) on the kernel 
of the operator Ai\i{u): 

m 

A,^{u)Au{xi, . . . , xjv; 2/1, . . . , yN-i) = A«(a;i, . . . , xjv; . . . , yN-i) (5.3) 

for m = 0,...,A^ — 1, {Ao{u) = 1 is implied). In the case m = and m = 1 this equality 
is obvious. For m = 2, . . . , — 1 it can be proved by induction using the relation (12. 7p . 
Indeed 

Am{u)Au{x; y) = iu + ihd^JAra_i{u)Auix; y) - e^™-^-^™A„_2(M)A«(a;; y) = 

m 

= {-{)--e^^^^'~''\u{x;y). (5.4) 
The formula (15. 3p for m = N — 1 and the relation (12. 6p give 

Cn{u)A^{x] y) = -e""^ AN^iiu)Auix; y) = 

= -{-tf-^e^-+^f-~^'^'^-y^^Au{x-y) = i-zf^'Au-^H{x;y), 

what means in turn the equality (ni). 

Note that if one sets m = in (15.40 then the term containing e^™"^™ does not arise 
and, consequently, we obtain zero in the right hand side. Thus we have proved directly 
that the kernel Au{x] y) satisfies to the equation (13. ip . 

Analogously, to prove (iv) we used 

Sjv(n) = iu-p^)B^_,{u) - e^^-'-^^BN_2iu), (5.5) 
and the action of operator Bm{u) on the kernel of the operator Ai\i{u) 

m k m 

Bm{u)A^{x-y) = e-^^i-ir~'Y.^-l)''+'l[ey^-'-^^ J] e^=-^=A„(x; y) (5.6) 

k=l 3=2 s=k+l 



17 



which can be also proved by induction. Exploring the relation fl5.5p and the formula (15.61) 
for m = — 1, — 2 we obtain 

N-l k N-1 

k=l j=2 s=k+l 

N-2 k N-2 

k=l 3=2 s=k+l 

N 

i=2 

Let us note that the properties (ni), (iv) up to the factors (— ^)^^~^ can be derived 
from (ll) using R-matrix formalism (see for example |KL3] . |KL2] . |KLlj ). Thes proper- 
ties lead to the Baxter equation for operator AAr_i(u). Indeed, taking into account the 
recurrent relations (12. 5p and Dn{u) = —c^'^Bn^u), property (ll) for u = v and applying 
the properties (ni), (iv) one yields 

iN{u)AN^i{u) = r^e-^^A7v-i(M - ih) + i^e'^'^AN-iiu + ih). (5.7) 



6 Conclusion 

The derivation of the eigenfunction of the open Toda chain can be considered as the 
fist step of the Separation of Variables problem for the periodic toda chain. As it was 
mentioned in the introduction the Separation of Variables is achieved by a special choice 
of the transition function. In the case of Toda chain this function should be chosen as 
follows [KLT] : 

?^£,7i,-..,7]V-i(^l? • • • '^^TV) = e ^T,,,...,T,^_^(xi, . . . ,X7V__i). (6.1) 

The eigenfunctions of the integrals of motion of the periodic Toda chain, which are the 
expansion coefficients for the operator t^iu), are represented in the form 

*b(xi, ...,xn)= J t/£;i,7i,...,7iv-i(3^i' • • • ,a;jv)<I>ij(7i, . . . ,'yN-i)lJ-il)d-f. 

IgiV-l 

The function (16. ip can be calculated by induction: having an expression for the (A^— 1)- 
particle transition function one can yield the expression for the A^-particle one. This 
method was proposed in the work |KL2] . The authors of |KL2j obtain the recurrent 
formula integrating over the variables 7^: 



/ JV iV — i \ 

(E Afc- E 7,>iv (6.2) 
^7i,...,7iv-i(a;i,---,a;jv-i)i^(A;7)/i(7)(/7, 
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where /i(7) is a Sklyanin measure described in section [Hand K{X] 7) is some kernel. This 
formula leads to the Mellin-Barns representation for the transition function. 

In the present paper it is shown that the recurrent integration can be realised in terms 
of the coordinates Xn (eq. fl3.20p ) . This leads in turn to the Gauss-Givental represen- 
tation (13.211) . which was obtained in |Giv] . [GKLO] from the other circle of ideas. The 
integration over the coordinates in the formula f l3.20p implies actually an action of the op- 
erator Ajv(7jv) on the A^-particle eigenfunction. The function (16. ip can be rewritten then 
in terms of the product of A-operators like (I3.16p . This results in the fact that the func- 
tion (16. ip inherits the properties of the A-operators discussed in section O For example, 
the Weyl-invariance of this function is encoded in the property (l). Due to the properties 
(ill) and (iv) the equation (15. 7p holds, and therefore the transition function (16. ip satisfy 
the Baxter equation 

where dk = {0, ... ,1, ... ,0) is a k-th base vector. Since the Sklyanin measure n{'y) has 
necessary translational properties the eigenvalue equation 

tiv(M)*i?(a;) = tNiu; ^)^ij(x), 

where 

N 
k=0 

is equivalent to the Baxter equation for the functions $e(7) 

t^hk, £^)$e(7) = + ^fiSk) + r^$,(7 - th6k). (6.3) 

This leads in turn to the occasion to present these function in the form 

N-l 

that is to the separation of variables (see details in |KLlj ). 

Availability of two kind of recurrent formulae - of type (16. 2p and of type (13.200 - is 
explained by the fact that the function ip.y{x) can be regarded as well as a function of x„ 
satisfying the differential equations (I2.16P and in other hand as a function of 7^ satisfying 
difference equations in 7^ [Bab], i.e. as a wave function of some dual model. The duality of 
the same kind appears in the Representation Theory. The infinite-dimensional Gelfand- 
Zetlin representation of Lie algebra qI{N) by shift operators in 7^- allows to obtain the 
Mellin-Barns integral representation |GKLj . while the Gauss representation of the same 
Lie algebra by differential operators in Xn leads to Gauss-Givental representation |GKLO] . 

We hope the method proposed for the XXX-model in |DKM] and developed here for 
the Toda chain (including the use of diagram technique) can be applied to other more 
complicated integrable systems. 
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